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Abstract 

We clarify questions related to the convergence of the OPE and conformal block 
decomposition in unitary Conformal Field Theories (for any number of spacetime 
dimensions). In particular, we explain why these expansions are convergent in a 
finite region. We also show that the convergence is exponentially fast, in the sense 
that the operators of dimension above A contribute to correlation functions at most 
exp(— a A). Here the constant a > depends on the positions of operator insertions 
and we compute it explicitly. 
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1 Introduction 

Conformal Field Theories (CFTs) appear in many branches of theoretical physics: fun- 
damental interactions beyond the Standard Model, string theory and quantum gravity, 
inflationary cosmology etc. They play a particularly important role in all problems involving 
renormalization group (RG) flows, such as the theory of critical phenomena. Indeed, a 
general Quantum Field Theory (QFT) can be thought of as an RG flow starting from a 
CFT in the UV and flowing to another CFT in the IR. 

CFTs are also under better theoretical control than general QFTs, with several tools 
which are not available or much less powerful without conformal symmetry. In this paper 
we will be concerned with one such tool: the Operator Product Expansion (OPE). 
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For simplicity we will be considering here the OPE of two identical scalar operators. 
In any Lorentz-invariant QFT, the operator product 0(x)0(O) in the short-distance limit 
x — > can be approximated by a sum of local operators with c-number coefficient functions 
depending only on x 2 : 

<t>{x)m X ~ ^(xV 1 • • .*«O M1 ... w (0) . (1.1) 
o 

The RHS will in general contain both scalars and operators of nonzero spin Z, and in the 
latter case their indices have to be contracted as shown. 

In the conformal case, Eq. (1.1) has three more specific properties. First of all, the x 2 
dependence of the coefficient functions is a power law fixed in terms of the operator scaling 
dimensions 1 

C(X 2 ) = Const. ( x 2)-A +(A o -/)/2 _ ( L2 ) 

Second, conformal symmetry allows us to classify all local operators of the theory into 
primaries, which transform homogeneously under the conformal group, 2 and their derivatives 
(descendants). We can now rewrite the OPE as a summation just over primaries, denoted 
by curly O: 

<t>{x)4>{y) = J2fnoP(x - y, d y )0(y) . (1.3) 
o 

The index contractions if O has nonzero spin are implicit. The coefficient function P is now 
a power series in d y and encodes the contribution of O and all of its descendants. Up to the 
overall numerical coefficient X^o > the form of this function is completely fixed by conformal 
symmetry in terms of the operator scaling dimensions. 

To make this more concrete, let us discuss the case of a scalar O in some detail. Consider 
the three point function (cfxfiO). On the one hand it is fixed by conformal symmetry: 

{<t>{x)ct>{y)0{z)) = f HO \x - y\ A °- 2A *\x - z\- A °\y - z\~ A ° . (1.4) 

On the other hand according to the OPE (1.3) it should be computable as: 

( ( j>(x) ( f>(y)0(z)) = U 4> oP(x-y,d y ){0(y)0(z)), (1.5) 

where 

(O(y)O(z)) = \y - z\~ 2A ° . (1.6) 

1 Strictly speaking, this is only guaranteed in a unitary theory so that the dilatation operator can be 
diagonalized. There is a special class of non-unitary theories called logarithmic CFTs in which the dilatation 
operator has Jordan blocks, and some correlation function contain logarithmic deviations from power-laws. 
In this paper we only consider unitary theories. 

2 = the group SO(D + 1, 1) of global conformal transformations of the Z?-dimensional Euclidean space. 
Our main results will apply for any D, although for D — 2 stronger results could be obtained using the full 
Virasoro algebra. What we call primaries is usually called quasi-primaries in 2D. 
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Demanding that the two expressions agree, the coefficient function P can be fixed unam- 
biguously. Here are the first few terms in the expansion: 



P(X,dy) = \X\ A ~ 2A * 



1 + ^ + ax li x v d^d v + (3x 2 d 2 + . . . , (1.7) 



2 

A + 2 A 

a = — — — , p 



8(A + 1)' 16(A-D/2 + l)(A + l) " 

For the OPE of two identical scalars considered here the coefficients depend on A = A^ 
but not on A^; in general they would also depend on A^ — A^ 2 . Such expansions have 
been worked out to all orders (and also for O of nonzero spin) in the 1970's [1-3] and show 
interesting structure visible already in (1.7). For example, the D dependence appears only 
in the terms multiplied by x 2 , and would therefore be sub leading on the light cone. Notice 
also that the D-dependent term becomes singular when Aq hits the scalar field unitarity 
bound D/2 — 1. This is not a problem since such an O is necessarily free and so f^o — 0. 

Once the OPE structure is determined, we can use it to express any n-point function as 
a sum of (n — l)-functions. Schematically: 

(<Kx)<Kv) II>(^)} = Y,f^P{x - y, d y ){0(y) J]>(^)> • (1-8) 

o 

For n = 3 there is a single exchanged primary O = ipi, and we go back to Eq. (1.5), but 
for n ^ 4 the sum will be infinite. Actually, it will be doubly infinite since P's are infinite 
series in d y . 

And here comes the third special property of the conformal OPE: it converges. By this 
we mean that the representations (1.8) are actually absolutely convergent at finite separation 
x — y, rather than being just asymptotic expansions in the limit x — > y. 

This property has two important applications: 

• Correlation functions of arbitrarily high order can be computed by applying the OPE 
recursively. Of course, to do this we must know all primary operator dimensions Aj 
and all OPE coefficients fijj. (collectively known as the CFT data). 

• Eq. (1.8) can also be used to constrain the CFT data itself, by means of an old 
idea known as the conformal bootstrap [4-6]. The point is that a conformal four 
point function can be computed using the OPE in three different channels: (12) (34), 
(13) (24), and (14) (23). That the results agree is a constraint on the CFT data. 

To explain how this works in the simplest setting, consider the correlator of four identical 
scalars 

{4>{x l )(f){x2)4>{x3)4>(x 4 )) = 2A ' 2A {Xij =Xi- Xj) , (1.9) 

X 12 X34 

constrained by the conformal symmetry to have this form with g(u,v) a function of the 
cross ratios 

2 2 2 2 

U=^ 5-, V = — 5-. (1.10) 

x 13 x 24 x 13 x 24 
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The same four point function can be computed by using the OPE twice in the (12) (34) 
channel, which gives: 

(0000) = T(fofP(xi - x 2 , d 2 )P(x 3 - z 4 , d A ) (0(x 2 )0(x 4 )) ee ^(/ ) 2 ■ ( LU ) 

O O X 12 X 34 

The conformal blocks Go{u, v) defined by the last equality (see [7] for a longer introduction) 
are completely fixed by conformal symmetry; they depend on O's dimension and spin as 
parameters. The function g(u,v) can thus be expressed as a series in conformal blocks: 

g(u,v) = l + J2(fo) 2 Go(u,v), (1.12) 
o 

where the first term represents the contribution of the unit operator in the OPE 0x0. 

Eq. (1.12) solves the problem of computing the four point function in terms of the CFT 
data. To derive this equation, we made an arbitrary choice of applying the OPE in the 
(12) (34) channel. However, the representation in the crossed channel (14) (23) should be 
equally valid. In the considered case of four identical scalar, the crossed and direct channels 
involve the same primary exchanges, and their equivalence can be expressed as a functional 
equation satisfied by the function g(u,v): 

g{u,v) = (u/v) A *g(v,u). (1.13) 

Substituting the conformal block expansion we get 

_ w a,j + J2(f ) 2 [v A *G (u 1 v) -u A *G (v,u)\ = 0. (1.14) 
o 

Since 1970's it was hoped that this equation and its generalizations to more general four 
point functions impose nontrivial constraints on the CFT data. In the last few years this 
hope is starting to be realized. A practical way of extracting information from Eq. (1.14) is 
as follows [8] (see [9-16] for generalizations, and [17] for a different method used for large 
N theories). 

Using conformal freedom, let's put the four operator insertions in a plane with three 
points fixed at 0, 1 and oo (see Fig. 1). Conformal blocks are then functions of the fourth 
point complex coordinate z, and the cross ratios are given by 

u = zz, v = (l-z){l-z). (1.15) 

Eq. (1.14) is then Taylor-expanded around the point z = 1/2, and each coefficient is equated 
to zero. This gives a system of linear equations for infinitely many unknown non-negative 
coefficients po = (fo) 2 - 

One can now try to make an assumption about the spectrum, like that there is no scalar 
operator below a certain dimension, or that there is a gap between the stress tensor and the 
next spin two operator, etc. Any such assumption can be encoded by saying that po should 
vanish in a certain interval. One then asks if there are nontrivial solutions satisfying the 
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Figure 1: Using conformal freedom, three operators can be fixed at x\ = 0, X3 = (1, 0, . . . , 0), 
Xi —> 00, while the fourth point X2 can be assumed to lie in the (12) plane. The variable z is 
then the complex coordinate of X2 in this plane, while z is its complex conjugate. 



assumption, which is a question of linear algebra and can be decided numerically. If there 
are none, then a CFT with such a spectrum cannot exist. 

For reasons which are not yet fully understood, this simple-minded approach gives pretty 
strong, sometimes optimal [9, 16], bounds on the operator spectrum and also on the OPE 
coefficients. Partly this success must be due to the fortunate choice of the Taylor-expansion 
point z = 1/2, halfway between the direct channel OPE limit z — > and the crossed 
channel limit z — > 1. But does the conformal block expansion converge near this point, as 
is necessary for the validity of the method? We noticed on several occasions that it does, 
supporting this statement with brief arguments of various convincing power. However, the 
issue is sufficiently important to warrant a detailed investigation. This is the main purpose 
of this paper. 

The discussion will be structured as follows. We begin in section 2 by sketching a well- 
known argument relating convergence of the OPE expansion to the convergence of a scalar 
product of two Hilbert space states in the radial quantization. We believe that this way of 
viewing the problem is correct, and will try to eliminate whatever doubts may still exist in 
its validity. First of all, to illustrate the Hilbert space techniques used in the argument we 
perform a few simple CFT computations and check that results agree with other methods 
(section 3). The heart of the paper is section 4, where we focus on a CFT four point 
function and discuss the rate (as opposed to just the region) of OPE convergence for this 
correlator. We show that this convergence is exponentially fast and derive the exponent. 
As a prerequisite to this result we derive asymptotics of the CFT spectral density weighted 
by the OPE coefficients, which is of independent interest. In section 5, we show that the 
convergence of the conformal block decomposition can be proven by the same methods, and 
that it is also exponentially fast. Consequences and potential applications of our results are 
discussed in section 6. 
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2 Radial quantization argument 



The following result concerning the convergence of the OPE is rather well known: expansion 
(1.8) will converge as long as x is closer to y than any other operator insertion: 

\x — y\ < min \zi — y\ . (2-1) 

i 

Sketch of proof . (see Polchinski [18], Sec 2.9) Let's quantize the theory radially with the 
point y as the origin. If the condition (2.1) is satisfied, we can find a sphere separating the 
points x,y from the rest of the operators (see Fig. 2). The correlation function in the LHS 
of (1.8) can then be viewed as the overlap 

(2.2) 

between the two states living on this sphere which are produced by acting with 0's and ip's 
on the radial quantization in and out vacua: 

|$) = <t>{x)<f>{y)\0) , (*| = (0| J] . (2.3) 

Furthermore, the state |$) can be expanded into a complete basis of energy eigenstates: 



•21 




'y i 

'Z3 




Figure 2: An n-point correlation function can be viewed as the overlap of two states in the 
radial quantization Hilbert space, living on a sphere. 



\Q) = Y,C n (x-y)\E n ). (2.4) 

n 

Since the radial quantization Hamiltonian is the dilatation generator, these states are 
generated by acting on the vacuum by local operators of definite scaling dimension A n = E n . 
Moreover, there is one-to-one correspondence between this expansion and the OPE (1.3). 
For every primary O in the OPE, the expansion (2.4) will contain a series of states produced 
by 0(y) and all of its descendants. Schematically: 

\E n ) = (d y ) n O(y)\0) E n = A + n. (2.5) 

The coefficients C n (x — y) with which such a state enters into the expansion is found by 
simply picking up the coefficient of {d y ) n in the conformal OPE. 
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Convergence of the OPE (1.8) now follows from a basic theorem about Hilbert spaces: 
the scalar product of two states converges when one of the two states is expanded into an 
orthonormal basis. Q.E.D. 

The above argument will provide a starting point for our discussion. Our goal will be to 
make it more explicit and quantitative, in particular by determining the rate of convergence. 

3 Operator formalism exemplified 

The purpose of this section is to provide some background material about the radial quanti- 
zation and the state-operator correspondence. This is pretty standard and may be skipped 
by the experts. 

3.1 Map to the cylinder 

One way to think about the radial quantization is by mapping the CFT from the Euclidean 
flat D- dimensional space to the cylinder K x S ^ 1 (Fig. 3). In the 2D case this is usually 
carried out by means of the logarithmic coordinate transformation. However, the map exists 
in any D because the cylinder is conformally flat: 

ds 2 cyl = dr 2 + dn 2 = r~ 2 (dr 2 + r 2 dn 2 ) =r' 2 ds\ D (r = logr , n 2 = 1) . (3.1) 

On the cylinder we have time translation invariance and the usual Hamiltonian quantization. 




Figure 3: The map between R and the cylinder. 



Going back to the flat space the cylinder time slicing is mapped onto the slicing by spheres, 
and we recover the radial quantization. 

It is a basic property of CFT that correlation functions on conformally flat backgrounds 
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can be computed by rescaling the flat space correlation functions. For scalars, we have: 3 

(<P(x) . . .)^ (x)dx2 = e-^^i^x) . . V > (3-2) 
which on the cylinder becomes: 

0(r, n) cyl = r A ^(x) R n . (3.3) 

Let us demonstrate this rule by computing the two point function on the cylinder starting 
from the flat space expression (1.6). We get [22]: 

(0(72, nOOtTi.nO)^ = e-^ll - 2 e ^( ni ■ n 2 ) + e- 2 ?\- A . (3.4) 

The cylinder times enter only via their difference = r 2 — ri, consistently with time 
translation invariance. It is also interesting to expand the two point function in the large f3 
limit: 

oo 

e-* A 5>e-*\ (3.5) 

n=0 

The powers of appearing in this expression are consistent with the tower of states (2.5) 
propagating along the cylinder. The coefficients c„ can be found explicitly: they are the 
Gegenbauer polynomials entering via their generating function: 

Cn = Cl*> (m-m). (3.6) 

We will reproduce this result via the operator methods below. 

3.2 Positivity and conjugation 

We next discuss reflection positivity and conjugations. On the cylinder we have the time 
reflection transformation 9 : r — > — t. Any (2ra)-point correlation function with operator 
inserted symmetrically under this transformation will be positive (in a unitary theory). In 
fact this correlation function can be viewed as computing the norm of a state generated 
by half of the operators (see Fig. 4) acting on the vacuum. This can also be seen in the 
Hamiltonian formulation where 0(r, n) cy i = e THcyl <p(0, n)e~ cyl , so that for an hermitean 
Minkowskian field 0(0, n) one has 

0(r, n)l yl = (e^0(O, n)e-^- 1 ) t = e'^^O, n)e T2f «* = 0(-r, n) cyl . (3.7) 

Reflection positivity in the Euclidean theory then corresponds to unitarity in the Minkowskian 
theory. For instance, {</>(— r, n)0(r, n)) cyl = (0(r, n)^0(r, n)) cyl corresponds to the (positive) 
norm of the state 0(r, n)|0) cy i. 

When we go back to the flat space, the transformation 9 maps to the inversion R : 
x — > x/x 2 . We conclude that a flat space correlator with operators inserted in an inversion- 
invariant way computes the norm of a state in radial quantization. Similarly, when dealing 

3 See e.g. [19, 20] for specific applications of this rule. The only exception from the rule is the stress 
tensor Weyl anomaly [21]. 
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Figure 4: Example of a reflection-positive six point function on the cylinder. 



with operators, hermitean conjugation corresponds to the inversion R. For real primary 
scalars, eqs. (3.3,3.7) then imply 

0( x )t = r - 2A *(j)(Rx) = R{(j){x)) . (3.8) 

(Analogous rules can be written for operators of nonzero spin.) For example, an inversion- 
invariant two point function can now be written as 

(O|0(x)V(x)|O), (3.9) 

where |0) and (0| are the radial quantization in and out vacua. 

Hermitean conjugation in the Euclidean theory extends to the conformal algebra gener- 
ators. In particular, we have 

Pl = K,. (3.10) 

This follows from the fact that a special conformal transformation is obtained from a 
translation by conjugating with the inversion: = RP^R. One way to check this result 
is by considering the corresponding differential operators expressed in terms of the cylinder 
coordinates: 

P M = -idp -> -ie -r [n M (9 r + (6^ - n^n. v )d/dn v \ , 

K,, = -i[x 2 d^ - 2x M (x<9)] ->■ -ie T [-n fl d T + (5^ - n^n^d / dn v ] . (3.11) 

We now see that they go one into the other under the time reflection. Moreover when going 
back to the Minkowskian cylinder, r — > it, we do find = PL implying the same relation 
for the Hilbert space operators. 4 

Relation (3.10) has several consequences. First, given a radial quantization state |$) 
and a conjugate state ($| the states 

P^|$> and (3.12) 

4 To avoid possible confusion, we should stress that the relation = P^ is only valid in the radial 
quantization, i.e. on the cylinder (Euclidean or Minkowski). For the original theory on flat Minkowski space 
the conformal generators are hermitean: IO = K, P^ = P. 
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will be also conjugate, as more general states obtained by acting with several P's and K's. 
Their scalar product will be positive as a norm should. 

The second application concerns expressing the operators <f)(x) and (f>{x)^ in terms of the 
radial quantization "asymptotic states", i.e. operators inserted at and oo. We have the 
usual QFT formula 

0( x ) = e ip - x (j)(Q)e- iPx </>{x)\0) = e ip ->(0)|0) = e iP ' x \(f>) , (3.13) 

where we introduced the notation for the state produced by acting with 0(0) on the in 
vacuum. On the other hand for <j){xy we have the following representation: 

0(ar)t = e iK - x 4>{oo)e- iK - x =^ (O|0(a;) t = (<P\e~ iK - x , (3.14) 

where the conjugate state (<f>\ is obtained by acting with 0(oo) on the out vacuum. The 
definition of 0(oo) as a conjugate of 0(0) contains as usual a rescaling factor necessary to 
get a finite limit: 0(oo) = lim^oo r 2A< *0(x). 

3.3 Reproducing two point function 

We would now like to illustrate the operator formalism with a few simple computations, 
checking that it gives results agreeing with other techniques. 

Let's begin with the scalar two point function. According to the above discussion, it can 
be evaluated in radial quantization as: 

(0(x 2 )0(xi)) = r^ 2A (0\O(Rx 2 yO(x 1 )\0) = r^ 2A {0\e- iK <^ 2) e iP ^\0) . (3.15) 

Expanding the exponentials, we get a sum of the same form as (3.5) with 

c„ = -^(^,n 2 |ri, ni > , \n,m) = (P-n 1 ) n \0), (n, n 2 | = (0\(K ■ n 2 ) n . (3.16) 

These are precisely the states mentioned in (2.5). The scalar products can be now evaluated 
using the conformal algebra commutation relations: 

[K(t, Pu] = 2t(5^D + M^u) , 
[M^, X p ) = -i{5 w X v - S vp X^ (X = P,K), (3.17) 
[D,P»\ = -iPp, [D,K tl \ = iK^ 

taking into account the action of the generators on the primary scalar state \0): 

K fl \O) = M^ u \O) = 0, D\0)=-iA\0). (3.18) 

The states \n, iij) with different n's are mutually orthogonal. 

It is not difficult to show that c n 's satisfy the Gegenbauer polynomial recursion relation 

nc n = 2(n + A - l)(m ■ n 2 )c n _! - (n + 2A - 2)c n _ 2 (3.19) 
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and therefore coincide with (3.6). For the first few terms the agreement is easy to check 
directly: 

c = (A|A) = l, ci = 2A( ni ■ n 2 ) , c 2 = 2A(A + l)(m • n 2 ) 2 - A , . . . (3.20) 

Notice that for ni = n 2 the states \n, ni) and (n, n 2 | become conjugate. The scalar product 
c n then computes their norm and must be positive. Eq. (3.20) is consistent with this. 



3.4 OPE and three point function 

A similar check can be done for the three point function. For simplicity we consider the 
case y = 0. Eq. (1.4) then predicts the following radial quantization matrix element: 

(Q\0{x 2 )(f>{x 1 )\(f>) = f^or2 2Ao r?°- 2A *\l - 2e~V ' n s) + e~ 2fi \-^l 2 . (3.21) 

If needed, this can be also expanded into powers of e _/3 = ri/r 2 with Gegenbauer polynomials 
as coefficients. 

We would like to see how the same matrix element is computed in the Hamiltonian 
language. The key is to interpret the OPE (1.3), (1.7) as an expansion of the state (p(x)\<p) 
into energy eigenstates on the cylinder, as we already remarked in section 2. Replacing 
d y —> —iP, we get 

0(*i)l0> = 5>^ r ?°" 2A * \ l + h P ■ n " r iX P ■ n ) 2 + ^> + --]\°)- (3-22) 
o 

For the bra part (0|(9(x 2 ) we use the same expression as in (3.15). 

The three point function computation is now reduced to evaluating matrix elements of 
products of K's and P's, just like in the two point function case. Using the values for a 
and (3 given in (1.7), it is easy to check explicitly that the first few coefficients of the e - ^ 
expansion agree with (3.21). 



4 Four point function 



In the previous section we saw how the state-operator correspondence can be used to 
compute two and three point functions. Since the form of these correlators is anyway 
fixed by conformal invariance, the point of this exercise was mostly to check the formalism. 
We will now apply these methods to the OPE convergence argument from section 2. 



4.1 Convergence 

Consider then the simplest higher point function: the correlator of four identical scalars. 
There are many ways to map this correlator to the cylinder, depending on the choice of the 
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radial quantization origin and the infinity. For simplicity we will first assume that two out 
of the four operators are inserted at these two special points: x\ = 0, X4 = 00. A more 
symmetric configuration will be useful later on. 

In this limit the four point function (1.9) maps to the radial quantization matrix element 

(0|0(x 3 )0(x 2 )|0>. (4.1) 

We wish to compute this matrix element as follows. First we expand the state 0(x 2 ) \<fi) by 
using the OPE in the 0(x 2 ) x 4>(xi) channel, as in Eq. (3.22). Second, by conjugating the 
same equation, we get an analogous expansion of the state (0|0(£ 3 ) into conjugate states 
generated by the action of fT's. Taking the scalar product, we get a sum of the form: 

00 

^(/^o) 2 ^e-^(0,n,n 3 |0,n,n 2 ), E n = A + n , . (4.2) 

O n=0 

Here \0,n, n 2 ) are states generated by acting on \0) with n powers of P with various 
contractions involving n 2 , while (O, n, 113] are the conjugate states involving K's and n 3 . 
For operators O which are Lorentz tensors these states also carry spin indices which are 
left implicit. The exact structure of these states is fixed by the functional form of the OPE 
coefficients, as demonstrated in Eq. (3.22) for the scalar O case, but we will not need to 
know it. 

The expansion parameter in (4.2) is related to the cylinder time interval: 

e -P = r 2 /r 3 & f3 = T3 -T 2 . (4.3) 

According to the criterion (2.1), the series is expected to converge whenever /3 is positive. 
This looks plausible as the terms with large dimensions are exponentially suppressed for 
> 0. However, what about the coefficients multiplying the exponentials? Can't they 
somehow overcome the suppression? Also, shouldn't one worry about the density of states? 

As we will now explain, convergence can be proven in a way which makes no assumptions 
about the size of the OPE coefficients or the density of states. To do this, consider the matrix 
element (4.1) as a scalar product of two states 

|$) = 0(x 2 )|0) and (tf I = (010(23) • (4.4) 

The norms of these states represent the correlators of the same form as (4.1) but at special, 
inversion symmetric, configurations. (If /3 > 0, then by dilatation invariance we can assume 
that r 2 < 1 < r 3 = r^ 1 so that the matrix elements defining the norms are radially ordered.) 
So these norms must be finite (and positive). In fact they can be computed by setting 
n 2 = n 3 in (4.2), which makes every single term positive, being a scalar product of two 
conjugate states. 

Now, a basic theorem about Hilbert spaces says that the scalar product of two finite 
norm states converges. The proof is by judicious application of the Cauchy inequality. We 
could apply it to the scalar product and bound it by the product of the norms. However, 
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this would not prove that the scalar product converges. For this we should consider the tail 
of the scalar product series (4.2) and show that it tends to zero as the cutoff E n ^ E* is 
sent to infinity. By Cauchy, such a tail is bounded by the product of the corresponding tails 
of the norm series. Since the latter tails consist of positive terms, they must tend to zero if 
the norms are to be finite. Q.E.D. 

With convergence of series (4.2) established, the next step is to study the convergence 
rate. 



4.2 Spectral density estimates 

It will be convenient to write the series as the Laplace transform of a "weighted spectral 
density" 

/>oo 

C(J3)= dEf(E)e- E ^ f(E) = y^p k 6(E-E k ). (4.5) 
Jo k 
Here the energies E/~ run over all the states present in the theory, primaries or not. The 
weights pk can be read off from (4.2). We will be considering the reflection positive case 
n 2 = n 3 for which all of these coefficients are positive. As explained above, the general case 
reduces to this one via the Cauchy inequality. 

For each fixed (3 > 0, the convergence rate of this integral is clearly controlled by how 
fast the spectral density grows at large E. We therefore begin by studying this second 
question. 

The key observation which will allow us to control the spectral density is that for /3 — > 
the four point function (4.1) can be computed by using the leading OPE for the product of 
two colliding operators <p{x^)(j){x2). Thus the small (3 behavior of the Laplace transform is 
fixed: 5 

C(f3) ~ /T 2A ^ {fi -)• 0) . (4.6) 

This asymptotics puts a strong constraint on the behavior of spectral density at large E. 
Were we to assume a power law growth of f{E), Eq. (4.6) would immediately fix the 
exponent and the prefactor: 

/(£) ~F(2A^ (naiVe) - (47) 

This of course cannot be interpreted literally as f(E) is a sum of 5-functions. But the 
integrated version of this equation is mathematically true: the integrated spectral density 



F(E)= I dE'f(E') (1.8) 

must behave asymptotically as 



o 



E 2A f 



5 Notation a ~ b means a/b — > 1 in the assumed limit. 
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This is known as the Hardy-Littlewood tauberian theorem and is surprisingly subtle to 
prove, if one wants to get the exact prefactor; see [23], Theorem 1.7.4 for integer-spaced 
spectra and Theorem 1.15.3 for the general case. The difficulty comes precisely from the 
fact that one does not make any regularity assumption about the weights pk apart from 
them being positive. 

Denote by Eul the energy above which the Hardy-Littlewood asymptotics (4.9) becomes 
effective, in the sense that the ratio F(E)/E 2A( i> can be bounded from above and below by 
positive constants. The precise value of Srl of course depends on these constants but we 
will not keep track of this dependence. Analogously suppose that the Laplace transform 
asymptotics (4.6) is effective for (3 ^ (3 . By dimensional reasoning, we expect E^ « l/(3 . 
An elementary proof of this relation can be given as follows. 



Start with the obvious estimate: 

£(/3) ^ / dE' f(E')e- E P ^ e~ E/3 F(E) . (4.10) 
Jo 

This is true for any (3. Fixing (3 = b/E with an arbitrary b = 0(1) we get 

F(E) ^ e b C{b/E) const. eV 2A <» E 2A * . (4.11) 

The second inequality is true as long as b/E ^ (3q, i.e. E ^ Ehl = A lower bound on F(E) 

can also be shown by similar tricks; see [24], Section 7.5.2. 

In the asymptotic limit E — > oo we could take const. — > 1 in the last inequality. The prefactor 
can then minimized by picking b = 2A^,, but it is still worse than the one in (4.9). To get the 
optimal prefactor, a much more elaborate argument is required. 



Finally, let us discuss the error term estimates in the Hardy-Littlewood asymptotics. 
The error in (4.6) is controlled by A , the dimension of the first non-trivial operator in the 
(j)X(f) OPE: 

C(j3) = /r 2A ^[i + o(/r Ao )]. (4.12) 

One could wonder if this power-suppressed error would translate into a similarly power- 
suppressed error term in (4.9). However, this expectation is wrong, proving once more that 
this theorem is more subtle than it seems. In fact, the best possible error estimate under 
these conditions is only logarithmic (independently of A ): 

E 2A <t> 

(4.12) =► F (E) = T{2A ^ + 1) [l + 0(l/logE)]. (4.13) 

See [23], Chapter VII, for this result and for the elementary examples showing that the 
estimate cannot be improved without further assumptions. 
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4.3 Comparison with partition function spectral densities 

To appreciate better the estimate (4.9), let us reproduce the parallel story for the usual, 
"unweighted", spectral density 

f (E) = J2^(E-E k ), (4.14) 

k 

which counts all the states with equal weight one. In this case the Laplace transform Cq((3) 
has the interpretation of the CFT partition function on S ' 1 times a circle of radius /3. The 
f3 — > asymptotics then takes quite a different form, as we expect 

C (P) ~ expt-vol^- 1 ) jr(/3)] (/3 ^ 0) , (4.15) 

where F{0) is the CFT free energy density in flat space IR D_1 x time. This is because in 
the high-temperature limit the curvature of S D ~ l can be neglected. 

Furthermore, by dimensional analysis the CFT free energy density must be given by (see 
e.g. [25]): 

= -K//3 - 1 . (4.16) 

The constant k > provides a measure of the number of degrees of freedom and one expects 
it to be O (central charge). In 2D, this can be shown rigorously using modular invariance: 

k = (tt/12)c (2D). (4.17) 



So, £q(/3) grows exponentially fast as j5 — > 0, and this requires an exponential growth of 
the corresponding spectral density. Using the stationary phase approximation, we find: 



f (E) « exp^ 1 



D 



7 



D 



vol(S 



D-l\ 



D - 1 



K 



I/O 



(4.18) 



This should be understood as an order of magnitude estimate after averaging; see [23], 
Section IV. 16, for a rigorous formulation. For D = 2 and k as in (4.17), this result reduces 
to the well-known formula for the asymptotic density of states in 2D CFTs. 

Eq. (4.18) puts our result (4.9) into nice perspective: the unweighted CFT spectral 
density entering the partition function grows exponentially fast, but the weighted ones 
relevant for the correlator computations have much slower power-like growth. This means 
that the OPE coefficients squared must be exponentially suppressed. 



4.4 Convergence rate 

We will now study the convergence rate of the OPE representation (4.2), written as the 
Laplace transform (4.5). In other words, we must study how fast the tail of this integral 

POO 

C{p,E«)= / dEf{E)e~ E r (4.19) 
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tends to zero for E* — > oo and a fixed > 0. 
The necessary estimation is done as follows: 



POO POO 

C(P, E.)=p dE [F(E) - F(E*)]e- E P ^ /? / dE F(E) 

JE t JE t 



e -E P 



Here in the first line we integrated by parts and then dropped the F(E*) term; keeping it is 
not useful because of the very weak error estimate in (4.13). We then used the asymptotics 
(4.9) and wrote the answer in terms of the incomplete gamma function. 

The most interesting case of this estimate occurs for E* 3> A^/ 0. Using the incomplete 
gamma function asymptotics, we get 6 

C(P, E.) < r(2A ^ + 1) £* 2A0 e-^ (E. » A^//3, E. » E HL ) . (4.21) 

We also added explicitly the condition that E* should be above Em,, so that the used 
spectral density asymptotics be applicable. 

Estimate (4.21) shows that the OPE expansion for the four point function converges 
exponentially fast. Given the importance of this result, we will give its alternative derivation 
directly from (4.6), bypassing the spectral density estimates. 



We can write: 

POO 

C(P,E*)= dEf(E)e- E ?' x e - E W-P), (4.22) 

J E, 

where f3' £ (0,/3) is a parameter to be fixed below. Now estimate the last factor by e~ E *^~^\ 
The remaining integral is C(/3',E*) but estimate it generously by £(/?')• We obtain 

Cifi, £*) < Ctfy-W-M < const. tfy*±+ e -E.V-P) f (4.23) 

as long as /3' is smaller than /3q introduced in section 4.2, so that the asymptotics (4.6) is effective. 
The /?' dependence is minimized for 

pf = 2A 4> /E, , (4.24) 

and we get 

C(P,E*)^CE* A *e- E *P (4.25) 

where C = const. (e/2A^,) 2 ^ is the same sub-optimal constant that we had in (4.11). For the 
estimate to be valid, we must have that j3' defined via (4.24) be smaller than both /3 and fio, which 
requires 

E* ^ 2A 4 ,/P, E* > 2A 4 ,/po • (4.26) 
These are the analogues of the two conditions present in Eq. (4.21). 

6 In this paper a < b means that limsupa/6 ^ 1 in the considered limit, E* — > oo in this particular case. 
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5 Conformal block decomposition 



In the previous section we discussed convergence of the OPE series representation of a 
four point function. This series (4.2) was a double series, the first summation being over 
primaries appearing in the OPE and the second over the descendants of each primary. As a 
matter of fact, this second sum over descendants computes the conformal block of a given 
primary: 

oo 

G (u, v) = J2 e~ E ^(0, n, n 3 \0, n, n 2 > . (5.1) 

n=0 

The values of u, v are found by taking the limit x\ — > 0, £4 — > 00 in their defining equations. 
We find: 

u = (r 2 /r 3 ) 2 , v = 1 + (r 2 /r 3 ) 2 - 2(r 2 /r 3 ) cos a , (5.2) 
where a is the angle between n 2 and n 3 . It is then easy to solve for z defined in (1.15): 

z = (r 2 /r 3 )e ta . (5.3) 

It's interesting to compare expression (5.1) with the standard definition of conformal 
blocks given in (1.11). These expressions are not manifestly identical, for the following 
reason. In (1.11) the usual OPE is used on both sides of the four point function. On the other 
hand, when deriving (5.1) we use the OPE on the RHS of the correlator, and the conjugate 
OPE on the LHS. One consequence of this difference is that expression (5.1) is term-by-term 
positive for real < z < 1, which corresponds to reflection positive point configurations, 
while the definition (1.11) is not manifestly positive. In spite of these apparent differences, 
the conformal blocks defined by (5.1) of course agree with the old definition as they should. 
As a check, we can evaluate the first few terms of small z expansion for the case of scalar 
O, using the OPE coefficients from (3.22) and norms from (3.20). Focussing on real z for 
simplicity, we find agreement with the formula 

GAfliz, z) = (^) ^ ,F 2 (f , f , ±=£ + 1; a±1 , A _ f + 1; (5 . 4) 

recently derived in [16]. 

Notice that representation (5.1) makes obvious two other general properties of the 
conformal blocks: they increase for real < z < 1 and decrease as z rotates away from the 
real axis. More precisely, the second property reads 

\G (z = re ia ,z = re~ ia )\ ^ G (r,r) (0 < r < 1) ; (5.5) 

it follows by Cauchy inequality. 

5.1 Convergence 

We now turn to the question posed in the introduction: what is the region and rate 
of convergence of the conformal block decomposition (1.12)? By convergence we mean 
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convergence of the partial sums 

lim V f 2 G (u,v), (5.6) 

A*— >oo z — » 

A ^A* 

and by convergence rate the speed at which the tails 

9a>a*( u > v )= fo G o(u,v) } (5.7) 

Aq^A, 

tend to zero. 

These questions can be answered based on the results of section 4. The point is that 
the tails of the conformal block series include but a subset of the terms in the tails (4.19) of 
the OPE series, identifying A* = E*. This means that the conformal block decomposition 
converges not worse than the OPE series. 7 Applying the main result (4.21) about the OPE 
convergence, we conclude that the conformal block decomposition converges exponentially 
fast for each z in the unit disc r — \z\ < 1 with the tail estimate: 

\9A^(z,z)\< AfV A * for A* > A^/(l - r) . (5.8) 



5.2 Extending the region of convergence 



Eq. (5.8) is sufficient to prove convergence of the conformal block decomposition at z — 
1/2, the point used in the conformal bootstrap applications described in the Introduction. 
However, this is not the best possible estimate of the convergence rate. This can already 
be suspected from the fact that conformal blocks are known to be real-analytic in the full 
z plane minus the cut (l,+oo) along the positive real z axis. We will now show that the 
conformal block decomposition converges in this larger region, with the following rate: 

|0a>a.(M)| < f74A7Tl) A * A ' |pW|A *- (5 ' 9) 



Here the function 



M = 7, , Z ,o (5-10) 



;i + \/r 



z 



maps the cut plane on the unit disc. Thus the convergence is exponentially fast everywhere in 
the cut plane. For \z\ < 1 we have \p(z)/z\ < 1 and so the rate of exponential convergence 
predicted by (5.9) is stronger than the one following from (5.8). For z = 1/2 we have 
p(z) w 0.17. 

As a matter of fact, the same function p(z) enters the individual conformal blocks 
asymptotics which can be derived [8] in D = 2 and D = 4 starting from the explicit 
formulas of [26-28]. While this by itself would not be enough to prove (5.9), it hints at the 
best possible character of this estimate. As we discuss below, this is indeed the case. 



7 To be precise, this argument should first be done for real z, where the series are term-by-term positive, 
and then extended to complex z by Cauchy inequality. 
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To show the new estimate, we will have to change the quantization origin. To compute 
conformal blocks at \z\ > 1 in radial quantization as an overlap of two states, we need to 
find a sphere separating x\ and x<i from x% (and X4 = 00). Such a sphere exists for all z in 
the cut plane, but for \z\ > 1 it can no longer be chosen centered at X\. 




Figure 5: Here we are considering four operators inserted at the shown points in a plane 
passing through the origin. We give complex coordinates in this plane. The circles are 
intersections of the spheres of radius 1 and r with the plane. The angle a is arbitrary. 



The needed geometry is best captured by inserting the four operator as in Fig. 5. On 
the cylinder, this configuration maps into the one shown in Fig. 6. Here ^34 are inserted 
at the cylinder time 0, while 01,2 at smaller time logr < 0. Computing the cross ratios and 
the value of z in this geometry, we find 

z = Ap/(l + pf (p = re ia ) p = p(z). (5.11) 

In other words, z(p) is precisely the inverse of the function (5.10). 

We now have to repeat the main steps leading to (5.8) in the new geometry. First of all 
we should express the radial quantization matrix element 

(O|0(x 3 )0(x 4 )</)(x 1 )0(a; 2 )|O) (5.12) 

using the OPE with respect to the new in and out vacuum positions. To begin with, the 
state 0(xi)0(x 2 )|O) is expanded into states generated by primaries and their descendants 
inserted at the origin. This expansion represents an equivalent form of OPE, the difference 
from (1.3) being that we have 0(^2 ) and not 0(y) in the RHS. It's trivial to work out 
coefficient functions in such modified OPE in terms of those appearing in (1.3). We will not 
need their exact form below. Schematically we will get: 

00 

<P( Xl )<P(x 2 )\0) = (2r)- 2A ^/^ ^r A - + "|0,n,a), (5.13) 

O n=0 

where the n sum is over descendants, the states \0, n, a) being generated by acting on \0) 
with n powers of P M in various contractions. Apart from the overall f^o factors, their form 
and the relative coefficients are completely fixed by conformal symmetry. In the considered 
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Figure 6: The configuration on S x R obtained from Fig. 5 by the Weyl transformation. 
The pairs of points ^3,4 and <f>\ t 2 are in antipodal positions on the spheres at the cylinder time 
and logr. Their positions on their respective spheres are rotated with respect to each other 
by angle a. 



geometry these states also have dependence on the angle a which we don't need to know 
explicitly. 

Conjugating the last equation, we can write an expansion of the state (0\<p(x 3 )<p(x^)\ in 
terms of the states (O, n, a — 0| obtained by acting with K's on the primary out state (0\. 
Taking the overlap of the two expansion, the four point function on the cylinder is expressed 
as the double series 

00 

J](/^o) 2 Yl rA ° +n (Oi ri, 0\O, n, a). (5.14) 

O n=0 

By the same argument as in section (4.1), we conclude that this series converges for all 
r < 1. 

To estimate the convergence rate we go to the reflection positive case a = and consider 
the limit r — > 1. In this limit the sum of the series (5.14) must behave consistently with the 
OPE. Two pairs of the operators are colliding simultaneously, so the asymptotics is given 
by ~ (1 — r)~ 4A< *. This in turn implies that the integrated spectral density must behave 
asymptotically as 

F(E)~—- 5.15) 

1 ; T(4A + 1) 1 ; 

By analogy with Eq. (4.21), the tails of the sum (5.14) obtained by cutting off the energies 
E n = A + n > E* must therefore decay exponentially as 

~ r(4A +i) E *** rB ' for ^» A ^/( 1 - r )- ( 5 - 16 ) 

This is the needed estimate, since for a = we have r = p(z) according to (5.11). We just 
have to interpret it for the conformal block decomposition. Just as before, conformal blocks 
are given by the sums over n in (5.14). Conformal block decomposition tails can be simply 
bounded by the OPE series tails with E* = A*, since the latter tails include more terms 
and all the terms in the considered case a = are positive. 
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Finally, consider nonzero a. The matrix elements entering the OPE series can be 
bounded by Cauchy inequality as follows: 

\(O,n,0\O,n,a)\ ^ ^(O,n,0\O,n,0)(O,n,a\O,n,a) = (O, n, 0\O, n, 0) , (5.17) 

where we used that the norms in the RHS are equal by rotation invariance on S D -\ In 
particular, individual conformal blocks can only decrease in absolute value as a becomes 
nonzero keeping r fixed: 

|CoU^o| ^ Go\a=o- (5.18) 

This inequality is the analogue of Eq. (5.5) in the new geometry. Therefore the conformal 
block decomposition tails at « ^ satisfy the same estimate as in the already considered 
case a = 0. This is precisely what is stated in (5.9), as r = \p(z)\ according to (5.11). 
Q.E.D. 

We now explain in which sense Eq. (5.9) is best possible. First of all, this can be checked 
numerically in simple solvable CFTs. For example, one can consider the four point function 
of free scalar boson in 4D whose conformal block decomposition is known explicitly [27]. 
Focussing on real < z < 1 for simplicity, one finds that the convergence rate is well fit by 

C(z)A:^p(z) A *. (5.19) 

Here j(z) ~ 1.5 weakly dependent on z. This shows that while the power-like prefactor 
is improvable, the exponential convergence rate is best possible. 

It is also instructive to see why it would be impossible to do better than p(z) A * using 
our method. Consider the following family of reflection-positive operator insertion configu- 
rations: 8 

x\ = r, Xi = 1, X2 = —kr, x 3 = —l/k. (5.20) 

For r < 1 and k between and 1 this family interpolates between the configurations 
considered in section 4 and here. We also have an interpolating family of four point functions 
expansions, with the operators in the RHS of the OPE always inserted at the origin and at 
infinity. The r — > 1 limit corresponds to colliding operators, and the Laplace transform has 
a power-like singularity. The order of this singularity will depend on whether k < 1 (one 
pair collides) of k — 1 (two pairs collide). But the leading exponential asymptotics of the 
Laplace transform tails at finite r will be the same: r A *. Now, the conformal invariant z 
for the configuration (5.20) is given by: 

z= {I + kfr / (I + krf . (5.21) 

To maximize the convergence rate for a given z, we must find a k such that the value of r 
found from this equation is smallest. It is then simple to check that k — 1 is always the 
optimal choice. 

8 A more general family of reflection-positive configurations obtained from (5.20) by rotating x\ and X4 
with respect to the origin by the same angle <p leads to the same conclusion. 
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6 Discussion 



The purpose of this paper was to provide a set of criteria for determining the region and 
the rate of convergence of the OPE and of the conformal block expansions routinely used in 
CFT. For simplicity, we focused on the four point function of identical scalars, although the 
method is general. Our main results are Eq. (4.21) and (5.16), establishing the exponential 
convergence rate of the OPE for two different schemes, depending whether <p{x)(j){y) is 
expressed as a sum of operators inserted at y or {x + y)/2. Convergence of the conformal 
block decomposition was shown to follow from that of the OPE. The optimal estimate 
(5.9) is obtained by using the second scheme. Near the point z = 1/2 used in the conformal 
bootstrap studies, truncating the conformal block decomposition at scaling dimension A* 
1 induces the error < 0.17 A *. 

Convergence is good news for the conformal bootstrap, because it means that the method 
and the results obtained so far are on mathematically solid ground. Fast convergence also 
means that approximate expressions for the conformal blocks obtained by throwing out 
the descendants beyond certain level may be used where all-order expressions are not yet 
easily available, like for external operators with spin. On the other hand, exponentially fast 
convergence also has a downside: it means that the very high-dimension operators effectively 
decouple, and it will be very difficult to learn about them by studying consistency conditions 
for the correlators of low-dimension operators, how it's done in conformal bootstrap. 9 

For CFTs with a gravitational dual, decoupling of high-dimension operators can be 
understood as an ordinary effective field theory decoupling of high-mass states in the bulk. 
In the AdS/CFT context, a CFT truncated to just the subsector of low-dimension operators 
has been dubbed "Effective CFT" [29]. Our results show that such a truncation gives an 
approximate description in any CFT, with or without AdS dual. 

To prove the exponentially fast convergence, we represented the four point function as 
the Laplace transform of a spectral density concentrated on the sequence of dimensions 
of the exchanged local operators (primaries and descendants). The weights were given by 
the OPE coefficients squared times kinematical factors (the descendant norms). The main 
observation was that the short-distance limit of the Laplace transform could be controlled 
by using the OPE for the colliding operators. By the Hardy-Littlewood tauberian theorem, 
this implied a power-like asymptotic behavior for the (averaged) spectral density. From here 
it was fairly simple to show that the Laplace transform tails at finite operator separation 
are exponentially small. 

The power-like asymptotics of the weighted spectral density is interesting in its own 
right. As we reviewed in section 4.2, the same density with unit weights (i.e. simply counting 
states) is expected to have an exponentially growing asymptotics ~ exp^i? 1- ^]. To make 
these two statements compatible, the individual OPE coefficients must be exponentially 
suppressed with the same exponent. 

Recently, /aster-t/ian-exponential suppression for OPE coefficients of high-dimension A 

9 One possibility to enhance the contribution of high-spin operators is to consider consistency conditions 
for the correlators analytically continued to the Minkowski space. 
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and not-too-high spin / operators was discussed in [30]. For large- N theories with gravity 
duals, they have an estimate 



This results follows from the fact that in the limit of large AdS radius the conformal block 
decomposition of the CFT four point function goes into the partial wave decomposition of 
the flat space gravitational S-matrix. High-dimension, low angular momentum contributions 
in (6.1) correspond to the transplanckian scattering at small impact parameters. Black 
hole production is then expected to dominate, while the elastic 2—7-2 amplitude will be 
exponentially suppressed by the black hole entropy. The corresponding piece of the four 
point function, f\ l times the conformal block, should then be similarly suppressed. 

Clearly, Eq. (6.1) comes from a very different physics input compared to our results. Thus 
it should not be surprising that it is stronger, predicting faster-than-exponential decay in 
the given range. 10 However, there remain points of contact between the two approaches. 
For example, it should be possible to show using our methods that, quite generally, the 
super-exponential suppression cannot extend for all spins (see below). This is also what is 
expected from the gravitational S-matrix point of view, since for large impact parameters 
near-elastic 2—7-2 scattering replaces black hole formation. 

Our results can be developed in several directions. Finding interesting generalizations 
involving different external states (unequal operators, operators with spin or with asymp- 
totically growing dimension) is one line of research. A less obvious question is how to make 
better use of the crossing symmetry constraint (1.14). As in [8], let us rewrite it dividing 
by the contribution of the unit operator: 



This "sum rule" has to be satisfied for any u, v where it converges (i.e. in the cut z plane 
as explained in section 5.2). A cornerstone of our argument was Eq. (4.6), obtained by 
using the OPE in the crossed channel. As such it is also a form of crossing symmetry, but a 
very weak one: it can be interpreted as saying that Eq. (6.2) must be satisfied in the limits 



However, Eq. (4.6) has a shortcoming: it does not say how fast the limit is reached. 
Suppose that we want to have C(/3) ^ 2/3 _2A< * for (3 < /3q, how small a f3o should we take? 
A priori this is not known, the answer may depend on a CFT and on a correlator under 
consideration. This ignorance feeds into the minimal value of the dimension cutoff above 
which our exponential asymptotic estimates become rigorous inequalities, as in Eq. (4.26). 
It would be nice to eliminate this non-constructive feature from our results, by showing that 
a universal /3q can be chosen. Preliminary work [31, 32] indicates that for D = 4 this may 
be done starting from the sum rule (6.2), at least for 0's of dimension less than about 1.7. 

10 For this reason it can be stated as a result for /a,; alone, since the conformal blocks behave like simple 
exponentials. 




(6.1) 




(6.2) 



o 



Z40,l. 
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The sum rule can be also used to explain why we believe (as mentioned above) that 
the super-exponential decay of all, rather than just low spin, OPE coefficients would be 
inconsistent with crossing symmetry. Morally such an assumption would be not too different 
from throwing out all high dimension primaries from the sum rule. Considering the case 
of real < z < 1, the functions Fq are peaked for z ~ 1 — 0(1/ A 2 ) and go to zero as z 
approaches 1 past this value. Thus it seems impossible to satisfy the sum rule near z = 1 
under such conditions. It would be interesting to work out the details of this argument. 

To conclude, we would like to emphasize once more the fundamental principle used in 
deriving the results of this paper: that the CFT correlation functions can be computed using 
the quantum-mechanical dynamics on the cylinder S ^ 1 x R, which is Weyl equivalent to the 
flat Euclidean space. It was assumed that we have a complete Hilbert space of normalizable 
states on the sphere S D ~ l , that the Hamiltonian can be diagonalized, and that the energy 
eigenstates form a basis. Given these assumptions it is possible to prove the existence of the 
OPE by expanding the state (f)(x)(f>(y)\0) into the states of fixed energy. It is also a simple 
matter to convince oneself that any state on the sphere gives rise to a local operator in the 
flat space (and vice versa), proving the state-operator correspondence. 

The discreteness of the spectrum per se was not assumed: all of our arguments go 
through if the sums over primary operators are replaced by integrals. The spectrum of 
operator dimensions is expected to be discrete in D ^ 3. In D = 2 there exist CFTs with 
a continuous spectrum of scaling dimensions, such as the free scalar boson with its vertex 
operators e ipX . The Liouville theory and non-compact sigma- models also have continuous 
spectrum. In such cases care must be taken in applying our results, since local operators 
generate states which are not normalizable (they are rather ^-function normalizable). 

As the very last comment, we should quote two monumental axiomatic quantum field 
theory papers [33, 34] devoted to the question of OPE convergence in CFT. Most of their 
efforts seem to go towards establishing rigorously the first step of our argument: that the 
OPE can be viewed as an expansion in the Hilbert space sense, after which the convergence 
is automatic. The question of convergence rate was not discussed. 

Acknowledgements 

We are grateful to John Cardy, Sheer El-Showk, Liam Fitzpatrick, Davide Gaiotto, Jared 
Kaplan, Zohar Komargodski, Hugh Osborn, Kyriakos Papadodimas and Leonardo Rastelli 
for useful comments and discussions. The work of S.R. is supported in part by the European 
Program Unification in the LHC Era, contract PITN-GA-2009-237920 (UNILHC), and 
by the Emergence-UPMC-2011 research program; he also thanks the Perimeter Institute 
for hospitality. The work of D.P. and R.R. is supported by the Swiss National Science 
Foundation under grant 200021-125237. 



25 



References 

[1] S. Ferrara, R. Gatto, and A. Grillo, "Conformal invariance on the light cone and 
canonical dimensions," Nucl.Phys. B34 (1971) 349-366. 

[2] S. Ferrara, A. F. Grillo, and R. Gatto, "Manifestly conformal covariant 
operator-product expansion," Lett. Nuovo Cim. 2S2 (1971) 1363-1369. 

[3] S. Ferrara, A. Grillo, and R. Gatto, "Manifestly conformal-covariant expansion on the 
light cone," Phys.Rev. D5 (1972) 3102-3108. 

[4] S. Ferrara, A. F. Grillo, and R. Gatto, "Tensor representations of conformal algebra 
and conformally covariant operator product expansion," Annals Phys. 76 (1973) 
161-188. 

[5] A. M. Polyakov, "Nonhamiltonian approach to conformal quantum field theory," Zh. 
Eksp. Teor. Fiz. 66 (1974) 23-42. 

[6] A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov, "Infinite conformal 
symmetry in two-dimensional quantum field theory," Nucl. Phys. B241 (1984) 
333-380. 

[7] M. S. Costa, J. Penedones, D. Poland, and S. Rychkov, "Spinning Conformal Blocks," 
JEEP 1111 (2011) 154, arXiv: 1109.6321 [hep-th] . 

[8] R. Rattazzi, V. S. Rychkov, E. Tonni, and A. Vichi, "Bounding scalar operator 
dimensions in 4D CFT," J 'HEP 0812 (2008) 031, arXiv : 0807 . 0004 [hep-th]. 

[9] V. S. Rychkov and A. Vichi, "Universal Constraints on Conformal Operator 
Dimensions," Phys. Rev. D80 (2009) 045006, arXiv : 0905 . 2211 [hep-th]. 

[10] F. Caracciolo and V. S. Rychkov, "Rigorous Limits on the Interaction Strength in 
Quantum Field Theory," Phys. Rev. D81 (2010) 085037, arXiv: 0912. 2726 
[hep-th] . 

[11] D. Poland and D. Simmons-Duffm, "Bounds on 4D Conformal and Superconformal 
Field Theories," JEEP 1105 (2011) 017, arXiv : 1009 . 2087 [hep-th]. 

[12] R. Rattazzi, S. Rychkov, and A. Vichi, "Central Charge Bounds in 4D Conformal 
Field Theory," Phys. Rev. D83 (2011) 046011, arXiv: 1009.2725 [hep-th]. 

[13] R. Rattazzi, S. Rychkov, and A. Vichi, "Bounds in 4D Conformal Field Theories with 
Global Symmetry," J. Phys. A44 (2011) 035402, arXiv: 1009.5985 [hep-th]. 

[14] A. Vichi, "Improved bounds for CFT's with global symmetries," JEEP 1201 (2012) 
162, arXiv: 1106.4037 [hep-th]. 

[15] D. Poland, D. Simmons-Duffm, and A. Vichi, "Carving Out the Space of 4D CFTs," 
arXiv: 1109.5176 [hep-th]. 



26 



[16] S. El-Showk, M. F. Paulos, D. Poland, S. Rychkov, D. Simmons-Duffin, and A. Vichi, 
"Solving the 3D Ising Model with the Conformal Bootstrap," arXiv: 1203 . 6064 
[hep-th] . 

[17] I. Heemskerk, J. Penedones, J. Polchinski, and J. Sully, "Holography from Conformal 
Field Theory," JEEP 0910 (2009) 079, arXiv: 0907. 0151 [hep-th]. 

[18] J. Polchinski, String theory. Vol. 1: An introduction to the bosonic string. Cambridge, 
UK: Univ. Pr., 1998. 402 p. 

[19] J. L. Cardy, "Is There a c Theorem in Four-Dimensions?," Phys.Lett. B215 (1988) 
749-752. 

[20] I. R. Klebanov, S. S. Pufu, S. Sachdev, and B. R. Safdi, "Entanglement Entropy of 
3-d Conformal Gauge Theories with Many Flavors," JEEP 1205 (2012) 036, 
arXiv: 1112.5342 [hep-th]. 

[21] A. Cappelli and A. Coste, "On The Stress Tensor Of Conformal Field Theories In 
Higher Dimensions," Nucl.Phys. B314 (1989) 707. 

[22] M. Luscher and G. Mack, "Global Conformal Invariance in Quantum Field Theory," 
Commun.Math.Phys. 41 (1975) 203-234. 

[23] J. Korevaar, Tauberian theory: a century of development. Springer, 2004. 497 p. 

[24] E. C. Titchmarsh, The theory of functions. Oxford, UK: Univ. Pr., 1939. 454 p. 

[25] S. El-Showk and K. Papadodimas, "Emergent Spacetime and Holographic CFTs," 
arXiv: 1101.4163 [hep-th]. 

[26] S. Ferrara, R. Gatto, and A. F. Grillo, "Properties of Partial Wave Amplitudes in 
Conformal Invariant Field Theories," Nuovo Cim. A26 (1975) 226. 

[27] F. Dolan and H. Osborn, "Conformal four point functions and the operator product 
expansion," Nucl.Phys. B599 (2001) 459-496, arXiv:hep-th/0011040 [hep-th]. 

[28] F. Dolan and H. Osborn, "Conformal partial waves and the operator product 

expansion," Nucl.Phys. B678 (2004) 491-507, arXiv:hep-th/0309180 [hep-th]. 

[29] A. L. Fitzpatrick, E. Katz, D. Poland, and D. Simmons-Duffin, "Effective Conformal 
Theory and the Flat-Space Limit of AdS," JEEP 1107 (2011) 023, arXiv : 1007 . 2412 
[hep-th] . 

[30] A. L. Fitzpatrick and J. Kaplan, "Analyticity and the Holographic S-Matrix," 
arXiv: 1111.6972 [hep-th]. 

[31] J. Espin, A study of the decoupling of high dimension operators in conformal field 
theories. Master Thesis, EPFL, 2011. 

[32] D. Pappadopulo and S. Rychkov, unpublished notes. 



27 



[33] M. Luscher, "Operator Product Expansions on the Vacuum in Conformal Quantum 
Field Theory in Two Space- Time Dimensions," Commun.Math.Phys. 50 (1976) 23. 

[34] G. Mack, "Convergence of Operator Product Expansions on the Vacuum in 

Conformal Invariant Quantum Field Theory," Commun.Math.Phys. 53 (1977) 155. 



28 



